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ABSTRACT 

Aims. We aim to study the geometry and kinematics of the disk around the Be star a Arae as a function of wavelength, especially across the 
Bry emission line. The main purpose of this paper is to answer the question about the nature of the disk rotation around Be stars. 
Methods. We use the VLTI/AMBER instrument operating in the K band which provides a gain by a factor 5 in spatial resolution compared 
to previous VLTI/MIDI observations. Moreover, it is possible to combine the high angular resolution provided with the (medium) spectral 
resolution of AMBER to study the kinematics of the inner part of the disk and to infer its rotation law. 

Results. We obtain for the first time the direct evidence that the disk is in keplerian rotation, answering a question that occurs since the 
discovery of the first Be star y Cas by father Secchi in 1866. We also present the global geometry of the disk showing that it is compatible 
with a thin disk + polar enhanced winds modeled with the SIMECA code. We found that the disk around a Arae is compatible with a dense 
equatorial matter confined in the central region whereas a polar wind is contributing along the rotational axis of the central star Between these 
two regions the density must be low enough to reproduce the large visibility modulus (small extension) obtained for two of the four VLTI 
baselines. Moreover, we obtain that a Arae is rotating very close to its critical rotation. This scenario is also compatible with the previous 
MIDI measurements. 



Key words. Techniques: high angular resolution - Techniques: interferometric - Stars: emission-line. Be - Stars: Keplerian rotation - Stars: 
individual (a Arae) - Stars: circumstellar matter 



1. Introduction 

The star a Arae (HD 158 427, HR6510, B3 Ve), one of the 
closest (d=74 pc, Hipparcos, Ferryman et al. ( 1997)) Be stars, 
was observed with the VLTI/MIDI instrument at 10 fim in 
June 2003 and its circumstellar environment was unresolved 
even with the 102m baseline (Chesneau et al. (2005), hereafter 
paper I), a Arae was a natural choice as first target due 
to its proximity but also its large mid-IR flux and its high 
infrared excess among other Be stars, e.g. E(V-L)~1.8 and 
E(V-12yL(m)~2.23. These first IR interferometric measurements 
indicated that the size of the circumstellar environment was 
smaller than predicted by Stee (2003) for the K band. The 
fact that a Arae remain unresolved, but at the same time had 
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strong Balmer emission, have put very strong constraints on 
the parameters of its circumstellar disk. Independently of the 
model, they have obtained an upper limit of the envelope 
size in the N band of (pmax- 4 mas, i.e. 14 R* if the star 
is at 74 pc according to Hipparcos parallax or 20 R* if the 
star is at 105 pc as suggested by the model presented in paper I. 

They finally propose a scenario where the circumstellar 
environment remains unresolved due to an outer truncation of 
the disc by an unseen companion. Nevertheless, this compan- 
ion would be too small and too far away to have any influence 
on the Be phenomenon itself. 

In order to study the inner part of this circumstellar 
truncated disk we have taken advantage of the higher spatial 
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resolution by observing at 2 //m with the VLTI/AMBER 
instrument in February 2005. It provides a gain by a factor 
5 in spatial resolution compared to VLTI/MIDI observations. 
We present in this paper these measurements showing, for the 
first time, a fully resolved circumstellar envelope in the Bry 
emission line and a clear signature of a Keplerian rotating disk 
around a Arae. We also discuss the challenging question on 
the nature of the geometry of the Be disks and particularly 
their opening angle since it is still an active debate. 

Following the Wind Comp ressed Disk model (WCD) by 
Bjorkman & Cassinelli lll993h . most authors have considered 
geometrically thin disks (half opening angle of 2-5 degrees) 
even if Owocki et al. (1996) have found that the equatorial 
wind compression effects are suppressed in any radiatively 
driven wind models for which the driving forces include a 
significant part from optically thick lines. Moreover, they 
found that gravity darkening effects can lead to a reduced mass 
loss, and thus a lower density in the equatorial regions. A 
wind compression effect is, however, not required to produce 
small opening angle of the disk. The investigation of accretion 
disks has shown that discs in hydrodynamical equilibrium and 
Keplerian rotation will not have much larger opening angles, 
since their scale height is governed by the vertical gas pressure 
only. For a disk to be thicker, either additional mechanisms 
have to be assumed, or it might not be in equilibrium at the 
radii in question (Bjorkman & Carciofi 2004). 

On the other side, Stee et al. (1995; 1998) claimed that Be 
disks must be more ellipsoidal in order to reproduce the strong 
IR excess observed and interpret the possibility for a Be star to 
change its spectral type from Be to B, and more rarely Be to 
Be-Shell type, i.e. where the disk is dense enough to produce a 
strong "shell" absorption. 

In the following we adopt, as a starting point, the same 
parameters for the modeling of a Arae, e.g. the central star 
and its circumstellar envelope used in paper I and summarized 
in Table 1. Following the polarization measurements Pl^O .6% 
and Position Angle (PA) of 172° by McLean&Clarke l 1 979l and 
Yudin hooH the disk major-axi s orientation is expected to be 
at about PA* 82° (Wood et al. ll996a[fl996b[ QuiiTenbach et 
al.U997). Assuming a stellar radius of 4.8 Rq and an effective 
temperature T^s = 18 000 K, the photospheric angular diame- 
ter is estimated to be 0.7 mas (Cohen et al. .1998.. Chauville et 
al.[2001 ). For the distance of 74 pc and a baseline of 60m at 
2 yum, Stee ( 120031) predicts the visibility of a Arae to be lower 
than 0.2, i.e. fully resolved. 

The paper is organized as follows. In Section |2l and |3l we 
present the interferometric AMBER observations and the data 
reduction. In Section |2 we try to obtain a first estimate of 
a Arae's envelope geometry using very simple "toy" models. 
Section|5]describes briefly the SIMECA code. In Section|6lwe 
present the best model we obtain with SIMECA that fits both 
the Bry line and the visibility modulus and phase as a func- 
tion of wavelength which allows us to infer the disk kinematics 
and its rotational velocity. Finally, SectionOdraws the conclu- 



parameter/result value 



Spectral type 


B3Ve 




18 000 K 


Mass 


9.6 Mq 


Radius 


4.8 Rq 


Luminosity 


5.8 IO'Lq 


Inclination angle i 


45° 


Photospheric density (p,,;,of) 


1.2 10"'^gcm-3 


Photospheric expansion velocity 


0.07 km s"' 


Equatorial rotation velocity 


300 km s-i 


Equatorial terminal velocity 


nOkms-' 


Polar terminal velocity 


2000 km S-' 


Polar mass flux 


1.7 lO^'M© year-' sr"' 


ml 


0.3 


m2 


0.45 


CI 


30 


Mass of the disk 


2.3 10-'°Mo 


Mass loss 


6.0 IO-'^Mq year-' 



Table 1. Model parameters for the a Arae central star and its 
circumstellar environment obtained from paper I 



sions from these first spectrally resolved interferometric mea- 
surements of a Be star at 2 jum. 



2. VLTI/AMBER observations 

Our observations of a Arae were acquired during a Science 
Demonstration Time (SDT) run, on the nights of 23 and 24 
February 2005, using the AMBER instrument on the VLTI 
in medium-resolution mode (R=1500). On the night of 24 
February, the observations were made with two UT telescopes, 
i.e. one interferometric baseline only, and consist of six expo- 
sure files, each containing 500 frames of 100 ms. On the fol- 
lowing night, three telescopes were used and three-baselines 
data were taken in a serie of three exposure files and another 
serie of two, each containing 500 frames of 70 ms. Without 
a fringe tracker, the integration time per frame must be short 
enough to minimize the smearing of the fringe visibility due to 
the beams jitter, while still having enough photons over the ele- 
mentary exposure. Its specific value is chosen depending on the 
atmospheric conditions. Immediately after observing a Arae, a 
nearby calibrator object was observed. 

3. Reduction of the interferometric Data 

The data has been reduced using the "ammyorick" package de- 
velopped by the AMBER consortium'. The principles of the 
AMBER data reduction have been described by Millour et al 
(2005). In addition to the tools furnished by the default pack- 
age, some specific treatment have been added for reaching an 
optimal precision on the interferometric observables. The vari- 
ous steps are the following. For each of the individual exposure 
frames, the complex visibilities are extracted from the various 
interferometric channels and calibrated using the photometric 

' see: http: //www-laog.obs. ujf-grenoble.fr/heberges/amber 
/article.php3?id_article=81 
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Date 


Timeobs 


Baselines 


Length 


P.A. 


Calibrator 






(UTC) 




(m) 


n 




(mas) 


24/12/2005 


08:46 


UT2-4 (Bo) 


80.9 


39 


HD124454 


1.52 ±0.02 






UT2-3 (Bi) 


46.4 


19 










UT3-4 (Bz) 


52.5 


81 






25/12/2005 


09:40 


UT2-4 (Bj) 


84.6 


52 


HD 124454 


1.52 ±0.02 



Table 2. Observation log, with the projected baseline lengths and angles. The diameter for the calibrator used is taken from the 
CHARM catalogue (Richichi & Pei-cheron l2002l) . 



channels and some internal calibration files. The piston p be- 
tween the beams is first estimated from the slope of the fringes 
with wavelength, and a correction phaser is applied on each 
frame. From this, the color-differential phase and visibility are 
calculated, for a given spectral channel A, with respect to a 
set of reference channel(s). We chose the reference channel to 
consist in the whole spectral bandwidth except the considered 
channel A. Due to the jitter of the beams and to the subsequent 
variations of output flux after the optical fibers of AMBER, 
not all the frames contain good quality fringes. A selection 
and weighing of the best frames above a given threshold is 
made based on a fringe quality criterion. This yields our color- 
differential estimators of phase and visibility with an optimal 
precision. For both the science source and its calibrator, the sta- 
tistical deviation of the differential phases over the frames serie 
is about 0.01 rad per spectral channel. The deviation of the nor- 
malized differential visibility is about 0.5%. These numbers are 
similar for all the baselines. 

The scaUng of the differential visibilities to "absolute" vis- 
ibilities is done through a best adjustment of the "science" data 
to the one measured on the calibrator star, whose theoretical 
visibility is known. For this, the histogram of the squared visi- 
biUties from all the science frames (with enough flux above the 
detector noise) is fitted, by a scaling factor, to the histogram 
of the calibrator star. With the data we have, this fit has been 
done globally, by integrating the spectral channels together for 
each frame. Therefore, some possible bias between the spectral 
channels would not be suppressed by the calibration process. 

In the present case, the observations of the nominal calibra- 
tor for a Arae were of poor quality, and another calibrator, at 
some distance on the sky (» 45°), had to be used. This does not 
seem to affect the calibration, though. The scaling fit between 
the histograms of the science and calibrator objects results in 
an error of about 5% on the visibilities, on every baseline. 



4. Study of the envelope morphology using simple 
models 

In this section we present the new AMBER data with the previ- 
ous MIDI data already presented in Paper I in order to obtain a 
first estimate of a Arae's envelope geometry using very simple 
"toy" models. 



S 10 



8 9 10 11 13 13 

Wavelength (inieron) 

Fig. 1. a Arae (unresolved star + uniform disc) model upper 
limit diameters (in mas), as a function of wavelength, from the 
2005 MIDI data described in paper I for the June, 16 (plain line) 
and June, 17 nights (dotted line). Note that the large values ob- 
tained for larger wavelengths are essentially due to a poor cali- 
bration and large error bars rather than a true physical effect. 



4.1. Envelope extension 

4.1 .1 . Extension in tine continuum at 2.1 jum 

Assuming that the visibility we measure in the continuum, Vc, 
is only due to the central star and its circumstellar disk we can 
write: 



Vc 



(1) 



where Vec and V^c are respectively the envelope visibility and 
flux in the continuum , V*c and F*c the star visibility and flux 
in the continuum and ¥,ot - Ffc+F^c. Since Vc is the measured 
visibility, in order to only estimate the envelope visibility, Vgc, 
we can rewrite Eq.^^ 



(2) 



The total flux is normalized, i.e. ¥ii,t-^ec+^*c-^- Since the 
star is almost unresolved 0.5<(/i*<0.7 mas which corresponds 
to 0.98< V*c<0.99 for the longest baseline at 2.1 jjm, we as- 
sume in the following that V*c=l. In order to estimate Vge we 
still have to determine the star and envelope contributions to 
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the 2.1 //m flux continuum. Since the envelope continuum flux 
in the visible and in the UV is negligible we have fitted the 
blue part of the SED using a T(.yy = 18000 K black body for the 
central star in order to deduce the envelope emission at larger 
wavelength. At 2.1/im we found that the star emission is still 
~ 1 .5 larger, i.e. 0.44 magnitude brighter, than the envelope con- 
tribution thus, following F*c=l-5 F^c, we obtain F*c=0.6 and 
Ffc=0.4. Equation|2lcan now be rewritten as: 

= 2.5K.- 1.5 (3) 

The continuum visibilities and envelope extensions we ob- 
tained following Eq. |3l assuming a uniform disk model, are 
given in table |3l 

4.1.2. Extension in tine Bry line 

We can define the same equation as Eq.|2for the envelope vis- 
ibility in the Bry line following: 



where V, and F,- are respectively the measured visibility and 
flux in the Bry line. and F,o/ are the quantities previously de- 
fined and Vf,- and Eg,, are the visibility and flux only due to the 
envelope, i.e. without the stellar contribution. We obtain using 
the Bry emission line profile F(,, =0.5 and F, =1.5 at the center 
of the Une. The values of V^,, we obtain using Eq.l^assuming a 
uniform disk model, are given in Table|3l 
Globally, we obtain that, using a simple uniform disk model 
for each measurement, the Bry emission in the line and the 
nearby continuum at 2.1 jim have the same extension about 7-8 
mas along the UT2-UT3 (Bi)and UT3-UT4 (B2) baselines and 
about 5 mas for UT2-UT4 (B3) and ~ 4 mas for UT2-UT4 (Bq) 
(which was measured at a slightly different PA). 

4.2. Comparison between MIDI and AMBER 
extensions 

In order to compare the extensions we obtain with AMBER at 
2.1 fj.m (see Tabled) with the MIDI data presented in paper 
I, we have calculated the envelope extension in the continuum 
between 7 and 13 jim following Eq.|2l assuming a uniform disk 
for each measurement. Unfortunately, these MIDI observations 
were conducted under unfavorable conditions with thin cirrus 
passing, and the data were reduced with one of the first versions 
of the software available. This led to eiTorbars between 8% and 
18%, insufficient to measure a significant change of the visi- 
bility and thus of the diameter between 7 and 13 ;um. Despite 
these limitations, these measurements have roughly given an 
upper limit value of the envelope extension around 10 mas in 
the N band (see Fig. 1) and are compatible with our AMBER 
measurements. 

Moreover, the extension of a Arae circumstellar environ- 
ment seems to be almost independent on the wavelength which 
puts strong constraints on the density law within the envelope. 
The large upper limit values obtained for larger wavelengths 
are essentially due to a poor calibration and large error bars 



rather than a true physical effect. For "classical" Be stars mod- 
els, where the density is slowly and regularly decreasing as a 
function of the distance from the central star, we predict an 
i ncre asing size as a function of wavelength as shown by Stee 
(1200 li) . A possibility to keep the same angular size between 2 
and 13 yum, may be a disk truncation by an external physical 
effect as proposed in paper 1. 

Using FEROS spectroscopic data and the fact that their 
VLTI/MIDl data are showing a nearly unresolved envelope 
they propose a possible disk truncation by a unseen compan- 
ion at a radius of about 154 R©, assuming a circular orbit for 
the companion with negligible mass. With 7?* = 4.8 Rq, this 
corresponds to about 32 stellar radii, i.e. 6.4 mas, which is in 
agreement with their estimate based on the VLTl/MlDl data 
for a disk truncated at 25 R*, i.e. 4 mas somewhat smaller 
than the companion orbit. This is also in agreement with what 
we have obtained with AMBER, i.e. an envelope extension at 
2. 1 fim around 7 mas which may appears larger than the 4 mas 
they obtained. Nevertheless, we must keep in mind that they 
estimate their extensions using a uniform disk for the star+disk 
emission whereas we use a unresolved star + a uniform disk 
only for the circumstellar envelope leading to a larger size. 

Finally, we have verified that our best model described in 
Section |5] was compatible with the MIDI data obtain in 2003 
showing a nearly unresolved envelope. 
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-10 1 , ^ , ^ , ^ ^ , ^ I ^ ^ I 

-10 -5 5 10 

(mas) 

Fig. 2. a Arae (unresolved star -1- uniform disc) model di- 
ameters (in mas) dependence on the baseline position angle. 
Full line: AMBER measurements in the continuum at 2.1 fim. 
Dashed lines: MIDI measurements. The length of the plot cor- 
responds to the values of the error bars (note the very large error 
bars for the MIDI data). The dotted line is the direction of the 
major-axis of the envelope estimated from polarization mea- 
surements (PA= 172°) obtained by McLean & Clarke (Il979l) 
and Yudin et al. (Il998h 
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Diameter = 7.6 mas 
Opening angle = 25° 
Inclination anlge = 4i° 









Disc : 

Diameter = 6.5 mas . 
Opening angle = 2° . 
Inclination anlge = 5: 







Envelope : 
. Diameter = 7 mas 
Opening angle = 120 
Inclination anlge 
Minor axis in the 
orientation of base 


■■■■■■■ 
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Disc : 

Diameter = 6.5 msS 
_ Opening angle = 2°' 
Inclination anlge = 
Polar jets : 
Diameter = 7 mas 
Opening angle = 50° 


7° 




Fig. 3. "Toy" models used in order to fit the AMBER and MIDI measurements. The 4 upper pictures are a cut of the circumstellar 
disk in a plane defined by the observer line of sight and the stellar rotational axis (the observer is on the right for each picture), 
the corresponding projections into the sky-plane with over-plotted the interferometric data points from MIDI and AMBER are 
the central pictures whereas a "3D artist view" is plotted into the lower raw for each model. 



Base n° 





1 


2 


3 


Description 


UT2-4 


UT2-3 


UT3-4 


UT2-4 


Length! (m) 


80.9 


46.4 


52.5 


84.6 


P.A. C) 


39 


19 


81 


52 


V, 


0.84 


0.80 


0.72 


0.73 


V, 


0.63 


0.70 


0.59 


0.55 


V 


0.60 


0.50 


0.30 


0.33 


v„. 


0.21 


0.50 


0.33 


0.19 


(^fc (mas) 


3.3 ±1.2 


6.8±1.6 


7.6±14 


4.5±0.9 


4>er (mas) 


5.3±2.0 


6.8±14 


7.3±2.1 


5.2±1.2 



Table 3. Visibilities measured in the continuum at 2.1 /im (Ve) 
and the Bry line (Vr) and deduced for the envelope contribution 
only, using a uniform disk model respectively in the continuum 
(Wee) and in the Bry line (V„). The corresponding angular di- 
ameters obtained for the envelope are given in the continuum 
(pec and in the line (per- 



4.3. Envelope geometry 

The envelope extensions presented in Table|3]are very sensitive 
to the sky-plane baseline orientation. This is particularly obvi- 
ous from Fig. 2 where we have plotted a Arae's (unresolved 
star -I- uniform disc) model diameters as a function of the sky- 
plane baseline position. In the following sections we present 
very simple models in order to constrain the geometry of the 
circumstellar envelope by fitting the data obtained (see Fig. 2). 
Note also that we present here models for the light distribu- 
tion on the sky, not "physical" models which will be presented 
hereafter. 



4.3.1 . Equatorial disk perpendicular to the polarization 

Our starting point is the equatorial plane position of the circum- 
stellar environment of a Arae deduced from the polarization 
measurement done by McLean et Clarke (1979) and Yudin et 
al. (1998). The polarization angle measured is RA.=172° and 
any flattened envelope model should have a semi-major axis 
perpendicular to this direction, i.e. around 82°. 
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In the following we consider simple axi-symmetric disk mod- 
els, presented in Fig.|3 with 3 free parameters: 

1 . The incUnation angle (0 between the observer and the polar 
axis (i=0° corresponds to pole-on). 

2. the opening angle (a) of the disk. 

3. the disk extension (a) in mas. 

The shape of the projection of the disk onto the sky-plane 
depends only on the two first free parameters, i.e. / and a. For 
all these "toy" models the observer is on the right for the 4 
upper pictures in Fig. |3 The coiTesponding projections into 
the sky-plane with over-plotted the interferometric data points 
from MIDI and AMBER are the central pictures whereas "3D 
artist views" are plotted into the lower row for each model. 

The ratio between the projected major-axis (a) and the 
minor-axis (b) of the envelope is given by: 



b cos / + Isirij sin (; - |) 

We assume that the major-axis of the envelope is oriented at 
P.A.= 82° (i.e. perpendicular to the polarization angle). Since 
B2 is oriented at 81° it is supposed to be a good estimate of 
the disk major-axis extension. In Fig. 0] we have plotted the 
ratio a/b as a function of the disk opening angle for different 
inclination angles between 35° and 60°. This angle range was 
determined following the inclination angle of 45° for a Arae 
determined from the fit of the circumstellar Wa, HyS and P/? 
emission lines in paper I. From Fig. 0] it is obvious that a 
precise determination of the inclination angle is mandatory in 
order to obtain an accurate opening angle estimation. 

Our first disk model 1 with /=45 °, a = 25 ° and a-1.1 mas 
corresponds to a projected ellipse with a a/b ratio of 1.2. The 
agreement with the observed data is good excepted for the B3 
and Bo baselines which present a smaller extension than pre- 
dicted. In order to fit the data for the B3 and Bo baselines we 
define a model 2 with i-51°, a -2° and 0=6. 5 mas, i.e. smaller 
and very thin compared to model 1 . In this latter case the a/b ra- 
tio is 1 .85 which allows to fit the B2 and B3 baselines and MIDI 
data but not the Bq and Bi baselines. As shown in Fig.|4]inodel 
(2), with a very thin disk, is not compatible with an inclination 
angle of 45° and we were obliged to use a larger inclination 
of 57°. From Fig. I^we can see that for an inclination angle 
of 45° the largest a/b ratio for a extremely thin disk (a ~ 0°) 
is only 1.41, i.e. V2. Moreover, changing the inclination an- 
gle will change the shape of the Ha, H/3 and P/3 line profiles 
which were used to infer the value of 45° obtained in paper 
I. Finally, using a simple disk model with various extensions, 
opening and inclination angles failed to fit simultaneously the 
AMBER and MIDI data especially for the Bo and B3 baselines 
which show a smaller extension along the PA = 39° and 52° 
sky-plane orientation. 

4.3.2. Polar-axis along the B3 baseline orientation 

Since the extension obtained along the Bo and B3 baselines are 
the smallest ones we may consider that these baselines are close 



2.0 



1.8 



1 .2 



.1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

Model 2 dnd 4 


















^1^45'\ 




^ '■'•1=40'' V \^^^ 


Model 1 


1=^5- ''■.. ^ ^ \^ 






20 40 60 80 100 

opening angle (degrees) 

Fig. 4. Flatness ratio of the projected envelope versus opening 
angle calculted with our simple disc model for different incli- 
nation angle. Horizontal lines show flatness from model 1 and 
model 2 

to the minor-axis of the envelope (i.e. along the stellar polar- 
axis). Thus, our model 3 with i=45 °, a = 120° and a=l.l mas 
is a very thick disk as shown in Fig.|4]with the polar axis close 
to Bo and B3 but thick enough to go through the Bi and B2 data 
points. This model is in fact very similar to the one presented 
in paper I but in this case it maybe diflicult to obtain a polar- 
ization large enough to be measured with a PA= 172° since the 
major-axis of this "disk" is not perpendicular to the polariza- 
tion direction. Nevertheless, the polarization was measured in 
the visible, and since it originates from the inner part of the 
disk, this model cannot be totaly excluded. 

4.3.3. Equatorial disk + polar enhanced winds 

One of the shortcomings of the previous models is that they 
cannot reproduce simultaneously the two main envelope 
characteristics : the polarization angle at PA=172° and the 
smallest extension along the PA= 39° whereas for the other 
AMBER baselines the disk is clearly more extended. 

Our last simple model 4 try to take into account these observa- 
tional characteristics by considering an equatorial disk oriented 
perpendicular to the polarization angle and flatened enougth 
to reproduce the diff'erence between the extension measured 
along the Bq and B2 directions. Moreover, in order to fit the 
extensions measured for the Bi baseline, we have added polar 
enhanced winds perpendicular to the disk as shown in Fig. |3l 
Thus, seen into the sky plane it is also possible to fit all the 
interferometric MIDI and AMBER data points. Nevertheless, 
due to the fact that the disk is geometrically very thin {a - 
2°) and for the same reasons already explain for our model 2 
the inclination angle must be around 57 °. Moreover, the po- 
lar enhanced winds must be very extended in latitude (with an 
opening angle about 50°) and dense enough in order to fit the 
Bi measurement. 
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4.4. Conclusion about the "toy story" 

Finally it seems that using very simple "toy" models we failed 
to reproduce the different extensions of the envelope as a func- 
tion of the baseline projection onto the sky-plane. Thus, in the 
following section we will use as a starting point our model 4, 
i.e. geometrically thin disk + polar enhanced winds which was 
the one with the better agreement with both AMBER data and 
polarization direction. Moreover, the formation of geometri- 
caly very thin disk, with a few degree for the opening angle 
seems to be the best scenario up to now, at least for the central 
Keplerian disk around Be stars. 

5. The SIMECA code: a brief description 

In order to constrain the physical parameters of the circumstel- 
lar environment of a Arae, we have used the SIMECA code. 
This code, described in previous pap ers (se e Stee & Arailjo 
Il994t Stee et al.[l995; Stee & Bittar l200ll) . has been devel- 
opped to model the environment of active hot stars. SIMECA 
computes line profiles. Spectral Energy Distributions (SEDs) 
and intensity maps, which can directly be compared to high 
angular resolution observations. The envelope is supposed to 
be axi-symmetric with respect to the rotational axis. No merid- 
ional circulation is allowed. We also assume that the physics of 
the polar regions is well rep resented by a CAK type stellar wind 
model (Castor et al. |l975l) and the solutions for all stellar lati- 
tudes are obtained by introducing a parameterized model which 
is constrained by the spectroscopic and interferometric data. 
The inner equatorial region is dominated by rotation, therefore 
being quasi Keplerian. The ionization-excitation equations are 
solved for an envelope modeled in a 410*90*71 cube. 

The populations of the atomic levels are strongly altered 
by non-LTE conditions from their LTE-values. For computa- 
tion, we start with the LTE populations for each level, and then 
compute the escape probability of each transition, obtaining 
up-dated populations. By using these populations as input val- 
ues for the next step, we iterate until convergence. The basic 
equati ons of the SIMECA code are given in detail by Stee et al. 
lll995l) . 

To take into account the photospheric absorption line, we 
assume the underlying star to be a normal B3 V star with 
Teff - 18 000K and R - 4.8Rq and synthesize the pho- 
tospheric line profiles using the SYNSPEC code by Hubeny 
(Hubenv .l98& Hubeny & Lanz 1995) . The resulting line pro- 
file is broadened by solid-body rotation and might be further 
altered by absorption in the part of the envelope in the line of 
sight towards the stellar disk. 

6. Using SIIVIECA for the modeling of a Arae 
circumstellar environement 

6.1. Fit of the SED 

Thanks to the previous study we have now a good starting point 
for the estimation of the global geometry of the disk around 
a Arae. Nevertheless, we want to use the SIMECA model al- 
ready described in paper I and in section 3 in order to obtain 
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Fig. 5. a Arae's Spectral Energy Distribution (SED) from var- 
ious sources in the literature (see text). Dashed line: emis- 
sion from the central star only assuming a black body with 
R*=4.8R0, T,//=18000K and d=105 pc. Fit of the SED with 
SIMECA taking into account the envelope free-free and free- 
bound contribution (plain line) between 0.3 and 20 fj.m 



a more physical scenario for this star. Thus, in order to ob- 
tain a SED which can be used to constrain our model we have 
collected photometric measurements from the largest spectral 
range available, i.e. fro m UV to mid-IR. We use d UV mea- 
surements from Jamar J 1,91^, Thompson Jl978l) . U magni- 
tude from Johnson ( 1966), BVRIJHK magnitudes from Ducati 
( 2002) . ISO data, MIDI SED measurements (2003) between 
7-13 fim and IRAS data. Of course, using non contemporary 
data for a variable star may be questionable but this is a good 
starting point especially for the shorter wavelengths related to 
the central star which is supposed to be stable, contrary to the 
circumstellar disk which may appear and vanish with a typi- 
cal time-scale of a decade and should modify drastically the IR 
excess of the SED. 

The SED we obtain from these various sources is plotted 
in Fig |5l The SED is dominated by the emission of the cen- 
tral star for wavelengths smaller than ljum which is assumed 
to be a black body with the following parameters from paper I: 
R*=4.8Ro, T^//=18000K and d=105 pc. We already discuss in 
paper I the fact that we were obliged to use a distance of 105 pc 
instead of the 74 pc obtained from Hipparcos measurements to 
fit the SED. Without considering any reddening, and keeping 
the Hipparcos distance, the radius of a Arae would be unreal- 
istically low (below 3.5 Ro) or the photosphere unrealistically 
cold (Tgff ~ 15000 K) thus we use the 105 pc determined from 
paper I in the following. The free-free and free-bound emis- 
sions from the envelope produce an IR-excess and dominate 
the SED for wavelengths larger than ~ 1 fim but the envelope 
remains optically thin. This emission depends mostly on the 
number of free electrons and on the temperature law in the cir- 
cumstellar envelope. As the envelope is almost fully ionized, 
we can consider that the global IR excess is only proportional 
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to the mass of the disk. In Fig |5]we present our best fit of the 
SED using a temperature law in the envelope with T(r) oc r"^^"* 
and a mass of the envelope of 4. 1 (± 0.5) 10""' Mq. 

The mass of the disk in SIMECA depends on six param- 
eters: three of them are related to the stellar mass loss: mass 
flux at the pole, Ci, and m^ The three other ones are related 
to the envelope kinematics: terminal velocity at the pole, at the 
equator and mi. 

The mass flux in SIMECA is given by: 



<I>(e) = <I)p„fe[l + (Cl-l)sin'«'(0)]. 



(6) 



where ml is the first free parameter which describes the varia- 
tion of the mass flux from the pole to the equator, and CI is the 
ratio between the equatorial and polar mass flux: 



CI = 



pole 



The expansion velocity field is given by: 



(7) 



10 20 30 40 50 60 70 

Latitude (decree) 



Fig. 6. Example of the SIMECA density distribution as a func- 
tion of the stellar latitude from the pole (9-0°) to the equator 
(6* = 90°) at 10 R,. 



vr{r,e) = Voiff) + [Voo(0) - Voiem - -y, 

r 



(8) 



We used y — 0.86 whic h is a typical value for early Be 
stars (Poe & Friend lll986l) : Araiijo & Freitas Pacheco C1984 : 
Owocki et al. lll994l) 1. 
with 



Voie) = 



l)sin'"H6')] 



PO 



(9) 



The second free parameter ml is introduced in the expression 
of the terminal velocity as a function of the stellar latitude: 



VUO) = Vcipole) + [Voo(eq) - Voo(pole)] sm"'^ie). 



(10) 




Fig. 7. Left: rotational velocities as a function of the stellar lat- 
itudes from the pole (0 - 0°) to the equator {0 = 90°) at various 
stellar radii (over-plotted on the graph). Right: terminal veloc- 
ity as a function of the stellar latitudes. 



The shape of the terminal velocity law as a function of the 
stellar latitude is plotted Fig.Qas well as rotational velocity 
for various stellar radii. 

Finally the density distribution in the envelope is given by the 
equation of mass conservation (see Fig.|6l: 



(11) 



The mass of the disk is obtained by integrating Eq. (11). It is 
possible to obtain a combination of these six parameters that fit 
the SED but the solution may of course not be unique. In order 
to discriminate between all the possible solutions we need to 
put more constraints on our modeling, namely try to fit the vis- 
ibility measurements in the lines and in the continuum for the 
various VLTI baseline orientations. 



6.2. Visibility modulus in the continuum 

In order to fit the visibility measurements described in section 
4.2 we simulate a thin disk + polar enhanced winds with 
SIMECA with a dense equatorial matter confined in the central 
region whereas a polar wind is contributing along the rotational 
axis of the central star. Between these two regions the density 
must be low enough to reproduce the large visibility modulus 
(small extension) obtained for the Bq and B3 VLTI baselines. 
Since p oc (see Eq. 11) we can build a model satisfying 
the above conditions by tuning both the mass flux and the 
expansion velocity in the envelope. The expansion velocity at 
the equator should be very small, i.e. a few km.s"', whereas 
at the pole it can reach larger values up to 500-2000 km.s"'. 
The density ratio between the equator and the pole must be 
around 10-100 to fit the data. The parameters we obtain for 
our best model are given in table l^and the 2.15 fim intensity 
map obtained in the continuum is plotted Fig |8] The largest 
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parameter/result value 



V sin i 


375 km S-' 


Inclination angle i 


55° 


Photospheric density (pphot) 


1.0 10"" g cm"'' 


Equatorial rotation velocity 


470 km S-' 


Equatorial terminal velocity 


1 km s ' 


Polar terminal velocity 


1000 km S-' 


Polar mass flux 


7 IO-'Mq year-' sr"' 


ml 


0.5 


m2 


100.0 


CI 


0.03 


Mass of the disk 


4.1 IO^'^Mq 


Mass loss 


1.3 10-*Mo year-' 



Table 4. Best model parameters for the a Arae central star and 
its circumstellar environment obtained from this work. Values 
for the other parameters that have not changed are listed Tabled 
and are from paper I. 



discrepancies between these parameters and those from paper 
I are that the star is now rotating close to its critical velocity 
(v=470 kms"', the disk is very thin but denser (by a factor 10)) 
and the expansion velocity is only 1 kms"'. The continuum 
map is seen with an inclination angle of 55°, the central bright 
region is the flux contribution from the thin equatorial disk 
whereas the smoother regions originate from the stellar wind. 
The brightness contrast between the disk and the wind is ~ 30 
but can reach 100 if you compare the inner region of the disk 
with the outer parts of the wind. This density contrast implies 
that we use a CI parameter smaller than 1, i.e. 0.03, in order to 
keep a suflicient density contrast equator/pole and a negligible 
equatorial expansion velocity (1 kms"'). 

The fit of the visibility in the continuum basically con- 
straints two parameters : the outer radius of the equatorial disk 
which will modify the visibilities for all the projected base- 
lines and the inclination angle which have a strong influence 
on the flattening of the projected equatorial disk into the sky- 
plane. The outer radius of the equatorial disk for our best model 
was set to 33 R* which is larger compared to the 25 R* found 
in Paper I but still in agreement with the 32 R* obtained for 
the companion orbit. Moreover, the inclination angle we ob- 
tain with our best model is 55+5° which is also larger but in 
agreement with the inclination angle of 45±5° found in paper 
I. 

6.3. Line profiles 

In order to put additional constraints on our modeling we also 
try to fit emission line profiles with SIMECA, namely Ha, H/3 
and Pa/?. These lines were also used in paper I but following 
a quite different scenario: a nearly spherical expanding and ro- 
tating envelope versus a thin disk + polar wind in this paper. 
Moreover, we also used in this paper the Bry emission line 
which is accessible to AMBER. We remind that the Ha and 
Hy6 line profiles were observed in April 1999 with the HEROS 
instrument at la Silla (Chile), the Pa^ profile was recorded in 
August 13, 2003 at the Observatorio do Pico dos Dias (Brasil) 




Fig. 8. Intensity map in the continuum at 2. 15 /im obtained with 
SIMECA for our best model parameters. The inclination angle 
is 55°, the central bright region is the flux contribution from 
the thin equatorial disk whereas the smoother regions originate 
from the stellar wind. The brightness contrast between the disk 
and the wind is globally ~ 30 but can reach 100 if you compare 
the inner region of the disk with the outer parts of the wind. 



and only the Bry line was contemporary to our AMBER in- 
terferometric observations. Knowing that a Arae is a variable 
star exhibiting line variations with typical time scales of a few 
months, we will only use the global shape of these line profiles 
in order to constrain the kinematics within the disk. We will not 
try to fit simultaneously the intensity of all these lines which is, 
moreover, not possible as shown in paper I. In our new model 
these emission lines originate from the dense and thin equato- 
rial disk and thus a fit of the shape of these lines will put strong 
constrains on the disk kinematics which is supposed to only 
rotate around the central star. Again, this is a quite different 
scenario compared to paper I where we used an expanding and 
rotating nearly spherical envelope to fit the lines. 

The observed Bry line profile is not showing a double- 
peaks structure due to the medium (1500) AMBER spectral 
resolution mode used for these first interferometric observation 
since the high (10000) AMBER spectral resolution mode was 
not available at this time. Nevertheless, we can see in Fig|9]that 
after a convolution with a 15 A gaussian, corresponding to the 
AMBER spectral resolution, the agreement with the modeled 
profile (lower row, right profile, dashed line) and the AMBER 
ones (plain line) is satisfying. For the other lines the fit of the 
peaks separation is in agreement with the keplerian rotation. As 
already discussed in paper I, we recall that between the 1999 
observations (Ha and Ufi profiles), the 2003 (Pa/?) and 2005 
(Bry) observations, we have decreased the density at the base 
of the stellar photosphere by 25 % in order to fit the line inten- 
sity. 
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Fig. 9. Line profiles modeled with the SIMECA code using 
a thin disk + polar wind scenario. Modeled profiles are dot- 
ted lines whereas observed ones are in plain line. For Bry we Fig. 10. Differential visibility across the Bry line profile for the 
have convoluted the modeled profile (dotted line) obtain with Bo, Bi, B2 and B3 baselines. The first picture from the top is 
SIMECA with a 15 A gaussian, corresponding to the AMBER the Bry line profile. The plain line are the fits we obtain with 
spectral resolution of 1500. The convoluted Bry profile we ob- SIMECA from our best model whereas the VLTI/AMBER data 
tain is the dashed line superimposed with the observed one are the points with error bars, 
(plain Une). 



Compared to paper I, the agreement between the modeled 
(dotted line) and observed (plain line) Ha line profile in Fig|5]is 
not as good. We are not able to reproduce the broad line wings 
which were mainly due to the nearly spherical expanding stellar 
wind used in our previous model. In our new scenario the geo- 
metrically thin and rotating equatorial disk produces a narrower 
Ha line profile. In order to obtain larger line widths we should 
have taken into account multiple diffusion that occurs preferen- 
tially in the line wings as shown by Poeckert & Marlborough 
ill979l) . This is especially true for Ha and Pa^S but less pro- 
nounced for Hj8 and Bry. Nevertheless, we concentrate our- 
selves to the double-peaks separation which is very sensitive 
to the rotational velocity law used and as already mentioned, 
the agreement we obtain was better using a keplerian rotating 
law within the disk. 
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6.4. Differential visibility modulus across the Bry 
emission line 

The differential visibilities curves are plotted Fig.^Jfor the Bo, 
Bi, B2 and B3 baselines. The agreement between the modeled 
visibility (plain-lines) and the VLTI/AMBER data is very good. 
Moreover, it was possible to reproduce the visibility decrease 
across the Bry line profile mainly due to the variation of the 
flux ratio between the unresolved star and the partially resolved 
circumstellar environment but also its shapes as a function of 
wavelength. For the B2 baseline mainly along the equatorial 



Fig. 11. Phase of the visibility across the Bry line profile for the 
Bo, Bl, B2 and B3 baselines. The plain line are the fits we obtain 
with SIMECA from our best model whereas the VLTI/AMBER 
data are the points with error bars. 



disk direction, the Bo and B3 baselines close to the equatorial 
disk, the curves, with a global "U" shape, present a visibility in- 
crease at the center of the line. This means that a Arae observed 
within a narrow spectral bandwidth of 15 A appears smaller at 
the center of the emission line compared to its extension in the 
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Fig. 13. y6 parameter as a function of the equatorial velocity tak- 
ing into account the stellar photosphere deformation due to its 
fast rotation. The plain line corresponds to physically possible 
cases whereas the dashed line corresponds to rotational val- 
ues larger than the "breakup" velocity. The dash-doted line also 
corresponds to unrealistic values but in this case we have not 
taken into account the stellar photosphere distortion larger than 
1.5 R* at the equator. Vm,„=351 kms"' corresponds to the min- 
imum vsini found in the lite rature for this star, i.e 288 kms 
by Uesugi & Fukuda ( Il970l) with i=55°. 



nearby emission line wings as seen in Fig^] On the contrary, 
for the Bi baseline oriented along the polar-axis the visibility 
curve present a single "V" shape decrease at the line center. 
These effects are mainly due to the kinematics within the disk 
and are confirmed by the shape of the phase of the visibility as 
we will see in the following section. 

6.5. Differential phase across the Bry emission line 

The shape of the differential phase across the Bry line profile is 
related to the photocenter displacement of the object as a func- 
tion of wavelength and is very sensitive to the /3 law used in the 
modeling as already shown by Stee ( 1996). For our best model 
we used a keplerian rotation within the disk and the resulting 
visibility phases as a function of wavelength for the Bq, Bi , B2 
and B3 baselines are plotted Fig. ^2 Again, the agreement be- 
tween the modeled phases and the VLTl/ AMBER data is very 
good which is the first direct evidence of tlie keplerian ro- 
tation witliin a Be circumstellar disk. Nevertheless, in order 
to know if this agreement is unique we have tested different 
models with different disk rotational velocities following: 

v^(r,0) oc sinel-j , (12) 

with j6 between the constant rotation (yS = 0) and the 
angular momentum conservation (J5 - 1 .0), respectively /3-Q, 
0.3, 0.4, 0.45, 0.5, 0.55, 0.7 and 1. Nevertheless, if you modify 
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Fig. 14. from our models and the AMBER data plotted as a 
function of the y6 parameter. xl,i„ is the minimum value obtained 
for a model with y6=0.48 which is very close to the Keplerian 
rotation (J3=0.5). 



the f3 law and still want to fit the double-peaks of the Bry 
emission line profile you need to simultaneously modify the 
equatorial rotational velocity of the star since the inner part of 
the disk, supposed to be in contact with the star's photosphere, 
rotates at the same velocity. Thus we were obliged to take into 
account the shape of the star photosphere and its distortion due 
to its fast rotational velocity as shown by Domiciano de Souza 
et al.(2003|) for the Be star Achernar which was rotating close 
to its critical velocity and thus exhibit a flattened photosphere 
with an equatorial vs polar radius ratio of 1.5. In Fig. I13lwe 
have plotted the "beta" law versus the equatorial rotational 
velocity that fit the double-peaks of the Bry line emission line 
profile, assuming a purely circumstellar rotating disk without 
expansion. This figure shows that if we assume that the stellar 
rotation is between v„„„ corresponding to the smaller measured 
vsini=288 kms"' by Uesugi & Fukuda (1970), with i=55°, 
i.e. v,„i„- 350 kms"' and v„,aj., corresponding to the critical 
velocity for this star, i.e. \,„„^=503 kms"'; the /3 parameter 
must lies between 0.33 and 0.54. This last value is in better 
agreement with the upper valu e of th e measured vsini=375 
kms"' by Bernacca & Perinotto 1 1910) with i=55°, i.e. v=457 
kms"' at the equator which may be a good indication that a 
Arae is rotating very close to its critical velocity. 

Moreover, in order to be more quantitative we have com- 
puted the total from our best model as a function of the p 
parameter. The results are plotted FigOand we clearly see 
that the better agreement {xl,j„) is obtained for a rotating law 
very close to the Keplerian rotation, i.e. xiiin~^-^^ ±0.04. 

6.6. Is the geometrically thin and rotating disk 
scenario unique ? 

Since our Keplerian disk -1- polar enhanced winds scenario 
seems to successfully reproduce all the available observables 
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Fig, 12. Intensity maps across the Bry line profiles within spectral channels of 15 A from which the differential visibility modulus 
and phase are estimated. Note that in order to increase the image contrast the central star and the continuum emission has been 
subtracted. 
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Fig. 15. a Arae differential visibilities (left) and phases (right) across the Bry line profiles for the Bq, Bi, B2 and B3 baselines. 
The plain line are the fits we obtain with SIMECA from a rotating and expanding scenario decribed in section l6^ whereas the 
VLTl/ AMBER data are the points with error bars 



(i.e. photometric, spectroscopic, polarimetric and interferomet- 
ric) we may wonder if this scenario is really unique, especially 
due to the number of parameters used in the SIMECA code. 
Thus, we have tested two other scenarios, the first one is based 
on the same global geometry of our best model, but we have 
now added an expansion component produced by an equatorial 
stellar wind. 

Since the peaks separation of the line profiles depends on 
both components of the velocity field, i.e. expansion and 
rotation, we have to decrease the rotational velocity used in 
our best model to keep the same observed peaks separation. 
Remembering that the lowest possible value of the rotational 
velocity is 351kms"' we can only use a maximum value of 
lOOkm-i for the expansion equatorial terminal velocity. The 



corresponding differential visibilities and phases we obtain are 
plotted Fig. 16.51 In this case, the decrease of the differential 
phase amplitude variations is due to the smaller value used 
for the rotating component. If the envelope was only in 
spherical expansion due to a spherical stellar wind, the phase 
variation across the Bry line would be zero since there is no 
photocenter displacement for a spherically symmetric velocity 
field. On the contrary, if you have a purely rotating disk, the 
photocenter displacement will follow the projection of the 
iso-velocities regions and will produce a typical "S-shape" 
as shown Fig. Thus, a decrease of the rotational velocity 
field component regarding the expansion component produces 
globally a smaller differential phase amplitude variation shown 
Fig. 14. The amplitude of the differential visibilty across the 
Une remains almost the same since it is less sensitive to the 
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Fig. 16. a Arae differential visibilities (left) and phases (right) across the Bry line profiles for the Bo, Bi, B2 and B3 baselines. 
The plain line are the fits we obtain with SIMECA from the model 3 decribed in section 4.3.2 and section l6^ whereas the 
VLTI/AMBER data are the points with error bars 



kinematics within the disk but rather to a global geometric 
change of the circumstellar environment and to a change in the 
star versus envelope flux ratio in the line and the continuum 
which remains unchanged by modifying the velocity fields. 
Nevertheless, the shape of the differential visibility is strongly 
modified and presents now a "V" shape as already mentioned 
in section 6.4. This shape variation is due to the fact that the 
expansion versus rotation ratio is increasing in the equatorial 
region and the difference between the pole versus equator 
kinematics is less pronounced when adding an equatorial 
expansion component. This analysis exhibit the importance of 
spectrally resolved interferometric measurements for the study 
of the kinematics within circumstellar disks. 



The second scenario is based on a quasi-spherical model as 
the model 3 already described in our "toy story" section. The 
kinematic used is very close to our thin disk + polar enhanced 
winds model, (i.e. Keplerian rotation with a 1km ' equatorial 
terminal velocity) but in this later case the disk is rather geo- 
metrically thick. Again, the corresponding differential visibil- 
ities and phases we obtain are plotted Fig. 15. Since the kine- 
matics remains mostly unchanged the fits of the differential 
phases are as good as for our best model but since the geometry 
of the disk is different the agreement with the differential visi- 
bilities is not as good. Again, thanks to the same disk kinemat- 
ics the shape of the differential visibilities is very similar but 
their amplitudes are not well reproduced. Nevertheless, these 
differences remain very small for an inclination angle of 55° 
and it is very hard to put an upper limit for the disk opening an- 
gle using this method, especially regarding the actual errors on 
the AMBER absolute calibrated visibilities. On the other side 
this study clearly shows that the equatorial region is a Keplerian 
rotating disk rather than an expanding wind. 



7. Discussion 

This study point out three important results touching lively de- 
bated issues: 



7.1. Keplerian rotation 

There were already some indications that the disk may fol- 
low the Keplerian rotation by other theoretical studies, for in- 
stance the results obtained by Hanuschik (il995i) regarding shell 
lines produced wi thin a K epleria n disk in hydrostatic equilib- 
rium. Hanuschik ( Il996l) : ( l200 Cf) also investigate the geomet- 
rical structure of the emitting part of circumstellar envelopes 
around Be stars and found a good agreement with spectro- 
scopic data using a thin disk in vertical hydrostatic and hori- 
zontal centrifuginal equilibrium, similar to a Keplerian accre- 
tion disk. These results were confirmed by Rivinius ( 1 19991) 
presenting high resolution echelle spectra for 6 B-type stars 
supposed to be seen edge-on and in good agreement with 
Hanuschick's models for the formation of shell lines in circum- 
stellar disks with Keplerian rotation. In a more recent paper 
Rivinius ( 1999) propose a scenario where the disk is no more a 
completely stationary structure but rather a succession of out- 
bursts which may form rings. But even within this scenario, 
part of the ejected matters attains sifficiently high angular mo- 
mentum to form a roughly Keplerian disk, at least for the im- 
mediate times of outbursts. Finally, another indirect argument 
in favor of a Keplerian disk is the success of the global oscilla- 
tion modeling already outlined in the paper review by Porter & 
Rivinius (2003). Thus, our results may be the way to put a final 
exclamation mark regarding the widely accepted fact that cir- 
cumstellar disk around Be stars are in Keplerian rotation. 
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7.2. Stellar rotation 

As already mentioned in section 6.5 we found that a Arae must 
be rotating very close to its critical velocity since we obtain 
— ~ 91%. This value is far above the conservative estimates 
of ~ 75 % usually found in the literature for Be stars. The fact 
that Be stars may be rotating much closer to their critical ve- 
locities than it is generally supposed was already outlined by 
Townsend et al. (2004) and Owocki ( 1996). This nearly critical 
rotation has quite profound implications for dynamical mod- 
els of Be disk formation and may be the clue for the Be phe- 
nomenon. It may bring sufficient energy to levitate material in 
a strong gravitational field or at least help other physical pro- 
cesses such as pulsation or gas pressure to provide sufficient 
energy and angular momentum to create a circumstellar disk. 
Moreover, observational evidences of this nearly critical rota- 
tion are growing such as the results obtained by Domiciano et 
al. (2003) using interferometric VLTl/VlNCl data of Achernar 
They measured a rotationally distorted photosphere with an ap- 
parent oblateness of 1 .56 which cannot be explained using the 
classical Roche approximation. This scenario follows the origi- 
nal picture by Struve of a critically rotating star, ejecting 
material from its equatorial regions. 

7.3. Polar wind enhancement 

Our interferometric measurements are evidencing a polar 
wind enhancement (see Fig. |8j which was akeady predicted 
for almost critically rotating stars. For instance, Cranmer & 
Owocki ( 1995) and Owocki & Gayley(J_99&) studied the ef- 
fects of limb darkening, gravity darkening and oblateness on 
the radiation driving mechanism and found that the tendency 
for the higher flux from the bright poles to drive material to- 
ward the darker equatorial region is outweighed by the oppo- 
site tendency for the oblateness of the stellar surface to direct 
the radiative flux to higher latitudes, i.e. a way fro m the equa- 
tor. The paper review by Porter & Rivinius (l2003h also outline 
the effect of the inclusion of nonradial line-driving force which 
reduces the effect of the wind compression to zero and, tak- 
ing into account the gravity darkening, results in a polar wind 
enhancement. This physical effect goes in the opposite direc- 
tion to the one proposed for the W ind Compressed Disk model 
from Bjorman & Cassi nelli |l993l) . In a recent paper, Kervella 
& Domiciano de Souza flOOff) have evidenced an enhanced po- 
lar wind for the Be star Achernar whereas this Be star presents 
no hydrogen lines in strong emission. Thus, it seems that a sig- 
nificant polar wind may be present even if the star is still in a 
normal B phase, i.e. this enhanced polar wind does not seem 
to be related to the existence of a dense equatorial en velope , as 
already outlined by Kervella & Domiciano de Souza ll2006l) . 

8. Conclusion 

1. Thanks to these first spectrally resolved interferometric 
measurements of a Be star at 2 /im we are able to propose 
a possible scenario for the Be star a Arae's circumstellar 
environment which consist in a thin disk + polar enhanced 
winds that is successfully modeled with the SIMECA code. 



2. We found that the disk around a Arae is compatible with 
a dense equatorial matter confined in the central region 
whereas a polar wind is contributing along the rotational 
axis of the central star Between these two regions the den- 
sity must be low enough to reproduce the large visibility 
modulus (small extension) obtained for two of the four 
VLTl baselines. This new scenario is also compatible with 
the previous MIDI measurements and the fact that the outer 
part of the disk may be truncated by a unseen companion at 
32 R*. 

3. We obtain for the first time the clear evidence that the disk 
is in Keplerian rotation, closing a debate that occurs since 
the discovery of the first Be star y Cas by father Secchi in 
1866. 

4. We found that that a Arae must be rotating very close to its 
critical velocity. 

5. These observations were done using the medium (1500) 
spectral resolution of the VLTl/ AMBER instrument and are 
very promising for the forthcoming AMBER high spectral 
resolution observational mode (10000) and the coupling of 
the Auxiliary Telescopes (ATs) on the VLTl array. 
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